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ABSTRACT 


Mays report deals with design optimization of Internal- 
Expanding Rim Brakes. A computer proaqram was developed to 
calculate the actuating force, torque, stopping time and 
drum temperature. fhe drum temperature is calculated by the 
finite difference method. 

A comparison of results has been made using a simplified 
equation that is in common use in engineering texts. 

Numerical optimization is shown to be a convenient tool 


for brake design. 
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SYYBILS AND ABBREVIATIONS 


PeeeeenGLISH LETTER SYMBOLS 


a Distance from pivot to the center of rotation (m). 
A Area of one lining shoe (m*). 

D Maden OL ELELCtION material (m). 

Be Biot modulus. 

Cc Specific heat (J/Kg-°C). 

Cc Boeemal capacity (37°C). 

d Distance from actuating force to the hinged pin (mn). 
Ge Rate of deceleration (m/sec? ). 

E Kinetic energy (J). 

f Frictional force (N) 

F Actuating force (N). 

Fy Fourier modulus. 

g Gravity constant (m/sect ). 

h Convection heat transfer coefficient (W/m? ~C) . 

k Thermal conductivity (W#/m-°C). 

M Friction moment (N-m). 

MH Norgnjal moment (N-m). 

N Normal force (N). 

p Pressure between lining and drum at any point (N/m*). 
BE Maximum pressure between lining and drum (N/m*). 
Q Heat generated (WW). 

c Inside drum radius (m). 

R Wheel radius (m). 

Bes. Thermal resistance (°C/¥W). 

t Tine (séc.) 

tk Thickness (m). 

fee Temperature (°C). 

I, Torque (N-m). 

v Velocity (m/sec.). 

VY, Volume (n°). 

W Vehicle weight (N). 





Be NOTATION 
Bay The thernal resistence between nede2e i and the 
adjoining node j. 


The temperatur2 of anode i at time step p. 


Ce muoee LETIER SY4BOLS 
9g The angle between the hinged pin and an element 
area on the lining. 
Q The angle at which the pressure between the lining 
and drum is maximum. 
Boeection coefficient. 
Cald triction coefficient. 
HOt rraiction coetficient. 
a Thermal diffusivity (m*/sec.). 


A Finite increment. 
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fen meOUU ST TON 


Brakes are mechanical devices for retarding the motion 
of a vshicle or machine by means of friction. Because of thea 
[itiarity of their functions, many clutches may also be 
ancluded here, assuming centrifugal forces are accounted 
EOL. 

A Simplified dynamic representation of a brake is shown 
Maeeig. 1. Two masses with inertias, I, and I, , rotating 
at the respective angular velocities w, and w, (one of which 
may be zero), are to be brought to the same speed by 
engaging the brake. 

The friction brake has three basic elements; two 
Opposing friction surfaces and a mechanism for forcing the 
friction surfaces into) contact. Whenever a friction brake is 
engaged *9 join two members having relative motion, there is 
a period of slip which may last several seconds. This slip 
is one of the chief merits of the friction brake; it absorbs 
siocks and prevents excessive torsicnal stresses on the 
power transmission system. On the other hand, siip is the 
iemeting £actor in friction clutch and brake performance; 
Betememeat iS generated in proportion to slip, torque 
transmitted, and period of slip. 

The following parameters ar= of interest in analyzing 
the performance of these devices; 

1. The actuating force. 

2. The torque transmitted. 
3. The temperature rise. 
4. The slip time. 

This report deals with Intarnal-Expanding Rin Brakes. 
This formulation also applies to internal-expanding clutches 


meecentrifugal forces are accounted for. 
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Ii, INTERNAL-EXPANDING 21M CLUTCHES AND BRAKES 


A. GENERAL MACHANICAL PRINCIPALS 

A brake or clutch assenbly, us2s a brake Shoe to which 
Memattach2d a friction material, called lining. The lining 
is riveted or bonded to the brak2 shoe as shown in Fig. 2. 
[The brake shoe is pivoted at a fixed point and the other end 
is subjected to a force which presses the shoe in contact 
with the irum. The force between the brake and the drum is 
radial as the drum rotates. If a point on the rotating drun 
surface first makes contact with the shoe at the end nearest 
the pivot, the shoe is termed a "trailing shoe". If it first 
makes contact at the other end the shoe is termed "leading 
shoe", the latter giving a higher braking torque than the 
former for a given braking force. 

The friction between the lining and the drum creates 
heat which is basically the conversion of energy of motion 
of the vehicle or machine to thermal energy at the friction 
Surfaces, namely the lining and the drum. This heat is then 
dissipatei and absorbed by the drum by conduction, 


convection and radiation into the atmosphere. 


B. FRICTION FUNDAMENTALS AND MATERIALS 

Friction mechanisms, such as brakes, are systems for 
converting mechanical energy into heat. Several basic 
factors affect friction and wear of materials used in brake 


Systems. fhe main factors are temperature, pressure, speed, 


Surface roughness, and type of anaterial. Some organic or 
molded friction materials show nos Changs in. fraction 
Characteristics with pressure, while others such as 


mempeered-metal materials dacrease in friction coefficient as 
pressure is increased. For metallic friction materials there 


is also a decrease in coefficient of friction as speed 


AZ 





increases. Tomperature effects upon the coefficient of 
friction vary widely with the type of materials used. 

In a tworshoe internal expanding brake there is a 
tendency for the brake drum to detorm under hard 
application. Drums odecome elliptical and the force to do 
mms 1S quite high and contributes to friction force. 

A brake or clutch friction material should have the 
following characteristics to a degree which is dependent 
upon the severity orf the service: 

mani gh and uniform coefficient of friction. 

2. The ability to withstand high temperatures, together 
with good heat conductivity. 

3. Properties which are not affected by environmental 
conditions such as noisture. 

4. Good resiliency. 

Beeeuiags resistance to wear, scoring and galling. 


ere BRAKE DRUMS 

One of the primary functions of a brake drum is that of 
absorbing and dissipating the heat developed during the 
application of the brake. A brake drum is a heat sink into 
which heat goes after it is created by the rubbing friction 
Meeone  Orake lining contact to drun. The brake shoe and 
lining permanently fixed on the axle, when actuated, 
contacts the drum under pressure to cause the friction to 
stop the vehicle. The enercy of motion of a vehicle is 
converted to thermal energy by the brake assemblies. A brake 
drum must have th2 capacity to absorb and dissipate this 
heat energy within the limits of the brake heat input. If 
this is not the case, the drum expands and the brakes fade 
or fail. The greater the mass of the drum, the more heat it 
Can absorb and store until such time as the heat can be 
dissipated by convection and radiation [Ref. 1]. 

An ideal brake drun would have the following 
@iaracteristics; 
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Merdtgn structural strength to resist bursting forces. 

pepeumators coefficient of friction. 

3. Hard surface to resist scoring. 

4. High heat conductivity *o rapidly conduct heat away 
from braking surfaces. 

5S. High emissivity factor to radiate heat from the drum 
surface to the atmosphere. 

6. High heat storag2 capacity to store heat from 
Successive brake applications until it can be 
dissipated. 

7. Good nachinability to permit boring of the drun. 


D. STATLC AND DYNAMICS ANALYSIS 
1. Assumptions 
In developing the equations, the following 
assumptions have been made; 

a. The pressure at any point on the shoe 1S proportional 
to the moment arm of this point from the pivot. 

b. The effect of centrifugal force may be neglected. 

c. The shoe is assumed to be rigid. 

Gaerne Eriction coefficient is a linear function of 
temperature and it does not vary with pressure, wear 
and environment. 

2. Pressure Concept 

Roranalyze an internal shoe refer to Fig. 2, which 
ShowS a shoe pivotei at a fixed point with the actuating 
force acting at the other end of the shoe. The mechanical 
acrangement does not permit pressure to be applied at the 
pivot, therefore the pressure at this point is zero. If the 
Shoe rotates through a small angle about A, the radial 
movement of any point on the arc of contact, is proportional 
£2 the moment arm of this point from the pivot. Assuming 
that the material of the pnrake lining and support obey 
Hooke's law, the pressure at this point will also be 
proportional to this moment arm. The distance is 
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Proportional to sin®d@. Therefore, the relations between 
pressure at any point and the maximum pressure, Bae will be 
given by the following formula; 
p - DP, 

sin6 sin6, (1) 
Foom this formula it can be seen that the frictional 
material at the heel, contributes very Little to the braking 
Se Lon , Eheretore it is bettec to begin the friction 
material at an angl2 6, greater than, say 0.15 rad. I+ can 
be seen also that the vressure will be maximum when 6 =90° 
or if the toe angle @, is less than 90°, then the pressure 
Will be naximum at the toe. For good performance it is 
recommended to concentrate as much frictional material as 
possible in the aesighborhood of the point of maximum 
pressure [Ref. 2]. 

See ActUuating Force and Torgue Caiculation 

From Fig. 2, it can be seen that the differential 

normal force on an element area of the lining will be; 


dN = pdA (2) 


whece dA is an area element of the lining and it's magnitude 
LS 3 

dA = rbdé CS 
In Equation 3, r is the inside drum radius and b is the drum 


Math, substituting for p and dA gives; 


p,brsing 


dN = ~ sing, dé (4) 


At the same point the dirferential frictional force is; 
df = udN (5) 
mete i is the coefficient of Friction. 
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mae actuating force, F, can be calculated using the fact 
that the summation of the moments about the hinge pin iis 
zero. The moment due to Frictional forces is; 


Bo 2 
M =9/ (r-cos8)df C6 } 


where a is the distance from the pivot to the center of 
Bocation. Substituting the value of dF and integrating 
eeom Gg, to 9, gives; 
2 
up, br 


aa = (en Oya Se a 
Me = sind, {(cos8, cos62)+5—(sin 8,-sin‘6,)} CZ) 


Bmere yy 1S assumed to be constant along the lining. 
Similarly the moment due to normal forces is given by; 


aoe 
MM = 
Me gf asin®8dN (8) 


Substituting the value Or aN and integrating 
from 0; to 92 gives; 
p,bra 
Mn = PUitmyenGace po ee cee a2) (9) 
The actuating force must balance the moments, therefore; 
MAM 
d Ci0) 





where dis the distance from the hinge to the point of 
Bomerecation of F. The torgue applied to the drum by the 
brake shoe is; 


To =,{°?rdf Gi 


Riter substituting the value of df and integrating ; 


up, br * 


To = (cos81-cos9> ) C2) 


sind 
a 
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4, Rate of Heat Generated andi Deceleration Calculation 

The differential rate of heat generated by an 

element area of the lining is equal to the velocity of the 

inside surface of the drum relative to the lining , times 

the differential frictional forse acting on the element 
ar ea; 


dQ = West (13) 


Assuming the brake is on a vehicle wheel with a radius of R, 
the inside surface velocity is egual to; 


vo — —V (14) 


where Vis the velocity of the vehicle and is a function of 
time. 

If vV=V(t) then Ne et) and the heat generated will be also 
a function of time. Suostituting the values of ve and df and 
integrating from @, to 6, , we get the following formula for 


the heat generated at any time t, 


Q(t) = Z)(cos81-coséz)V(t) (15)) 


sind, ( 
The kinetic energy of a vehicle of weight W is given by; 


i W 
E = me) C16 ) 


Note that if the brake is ona four wheel vehicle, ther 


o 


Will be eight shoes. Assuming all are leading shoes, each 
Will stop one-eight of the vehicle weight, so W/8 must be 
used in Equation (15). The rat2 of change in the kinetic 
energy is; 

dE 


ac = (SY aE (17) 


From the energy conservation law the rate of change in the 
kinetic energy is equal to the heat generated; 


te 





Q(t) = = (18) 


p@pstituting the value of Q(t) and dE/dt, it is seen that 
the velocity V(t) cancels and so the deceleration is nota 
function of time. Therefore the deceleration, dc, is; 


p,bu 
sind 
a 


de = we (=) 


2 
ae (=) (cos@1-cos8 2) (1) 


The velocity at any time is; 
V = V,-det (20)) 


where V, is the initial velocity. Substituting the velocity 

il Equation (16), yields the rate of heat generated asa 

honction of time, 

p,bu 

sing 
a 


Q(t) = (Z)(cos6 -cos82)(V,-det) C21) 


In this study the friction coefficient was taken as constant 





up toa temperature of 90°C and after 90°C, decreases 
linearly to zero at a specified temperature, er 
Uo Ree e 
Hoh 
w= 4 Ur pp (T-90) OR <a (22) 
0 one Tmax. 


where Ho iso the cold Coeritcient Of EYiction and Un is the 


moe coefficient of friction. 


E. SURFACE TEMPERATURE CALCULATION 

Samce the function of a brake is to convert kinetic 
energy into heat, surface temperatures of brake linings and 
dcums are nost important. Therefore it is necessary to know 
the temperature of the mechanism during and after any stop. 
The temperatures were calculated by the finite difference 
method. 
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1. Assumptions 


a. One dimensional heat flow-The heat flow is from the 
inner surface to the outer surface of the druna. 

b. Constant heat transfer coefficient. 

c. No heat dissipated by radiation. 

d. The heat is generated on the inner surface. 

Eee Lemperature Analysis 
ae | 6TNeOLY 
The differential equation to be solved in order 

to find the temperature in the drun, based on the 
assumptions, is; 


fee 2. 1, aT 23 
ox? * & (>? t 2 


with the following boundary conditions: 

at x=0 heat is generated, 

at x=tk heat is transfer2d to the atmosphere by 

convection. 
In the eguation above k is the thermal conductivity, a is 
the thermal diffusivity, tis time and tk is the drun 
thickness. This eyuation can be solved by the finite 
difference method "Ref. 3]. The finite difference model 
used here is shown in Fig. 3. The rate of change with time 


of the internal energy of a node 1 is approximated by; 
ptl_»p 
AE i A 


= (24) 


AG 


where p is the density, cis the specific heat and V, is 
the drum volume. 


Now define the thermal capacity as 


= 20 
Cc, p,c,AVo, (25) 


The forward difference equation for all nodes and boundary 


Sonaitions is; 
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Te 7P pbt1_pp 


p i at a6 
ae +l = C.—— (26) 
Reh. ij 1 AG 
where is the thermal resistance 
“th, ij +1 ; 
Solving the above eguaticn for a gives; 
pP 
i. ; 
ree? = (qP+ Jd orc neon ae (27) 
ar foo. Ly 


The thermal resistanc? can be calculated from the geometry 
and boundary conditions (Ref. 3]. To ensure 
Stability At must be equal or less than the following 
nodal relation; 





Cy 
aa < can i (28) 
Ren ig 
With the assumptions made, the drum can be viewed as an 


infinite plate, with heat generated at the surface of the 
first noise, as shown in Fig. 3. It is assumed that in every 
drum, there are two shoes and that both ar2 leading shoes. 
Therefore, two times = must be taken. 
pP 
+ 
Te 


i = (2Qh+ = _—e— nie (29) 
igh — et ci CR meh 


Se FOcMulation 

In the computer program 5 nodes were taken. In 
Bea@er tO sheck accuracy, the program was run with 7 and 10 
nodes. In each case the result was the same within 5 °c. 
The heat is generated in the inner drum surface. Therefore Q 
appears in the formula of temperature in the first node and 
for all the other nodes Q is equal zero. With the 
assumptions mentioneli above, the heat transfer through the 
drum is solved as a heat “transfer problem through an 
infinite plate, with heat generation at the inner surface 
and with a heat convection boundary on the outer surface as 
Srown in Fig. 3. 2quation (29) can be simplified using two 


dimensionless parameters, Biot ani Fourier modulii, 
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hAx 
aL k 


wD 
iI 


(30) 


mar . aAt Cel) 
(ax) 

The final equations for calculating the temperatures at the 

nodes now become; 


For the first node; 


2QP at 
prt ah Z p p 
eS eee ral © 27 0 (32) 


For the interior nodes; 


eee Pp Pp Ds p 
fee Pott, y+ Ge - 2T;} C32) 
For the Last node; 
pat _ p Lt ep _1)\mP 
Be ~ CPotT, ,t8,1.+ (gpy -By-1)T,) (34) 


Es BRAKE DUTY CYCLE 

In addition to tn2 parameters mentioned above the design 
of a brake depends on the initial speed, final speed, number 
of stops, and the rest time between each stop. In this 
analysis a general duty cycle was considered so that the 
initial speed, final speed ani the acceleration period 
between stops can be different for each part ot the design. 

In the design examples presented here, a vehicle was 


stopped four consecutive times with the following cycle; 


Initial Final Rest 

Speed Speed 

m/sec. n/sec. sec. 
1 Zo. 0 0.0 20.0 
2 25. 0 Q.0 20.0 
3 ZO O70 207310 
4 Zo) 0.0 = 





eter lui aaAl TON 


A. INTRODUCTION 

Engineering analysis using the digital computer has 
become commonplace. It 1s less common to use the? computer 
t> make the actual design decisions, such as sizing of 
structural members or placement of mechanical linkages. This 
may be largely attributed to the fact that fully automated 
design reguires techniques that are unfamiliar to much of 
the engineering community. 

ioe any engineering problems, it is necessary to 
d2termine the minimum or maximem of a function of several 
variables, linited by various linear and nonlinear 
inequality constraints. It is seldom possible, in practical 
applications, to solve these problems directly, and 
iterative methods are used to obtain the numerical solution. 
Meacawine calculation of this SOUL Ont, of ‘colizse, 
d2sirable. The CONMIN program is available to solve a wide 
variety of such oroblems { Ref. 4}. 

CONMIN is a FORTRAN program, in subroutine forn, for the 
minimization of a amulti-variable function subject to a set 
of inequality constraints. The basic optimization algorithn 
is the Method of Feasible Directions [Ref. 5]. The user must 
provide? a main calling program and an external routine to 
EeMeart@ate the objective and constraint functions and to 
provide gradient information. If analytic gradients of the 
objective or constraint functions are not available, this 
information is calculated by finite difference. While the 
Meogram iS intended primarily for efficient solution of 
constrained problems, unconstrained function minimization 
problems nay also be solved, and the Conjugate Direction 
Method of Fletcher and Reeves is used for this purpose 
Peet. 6}. 





B. DETLALTION OF TERMS 

Most lisciplines have a unique set of nomenclature used 
to describe the concepts within that discipline. Sone of the 
commonly used terms in numerical optimization are summarized 
here. 

Jbjective- The value of the function which is to be 
Minimized or maximized during the optimization process. 
Synonyms are cost, merit and payoff. The common mathematical 
designation is F(X). In the present study the objective was 
to minimize the material in the brake drun. 

Design variabl2s- Th2 parameters to ode changed 
during the optimization process in order to minimize or 
Maximize the value of the objective function. Synonyn; 
decision variables. The common nathematical designation is 
the vector xX. Design variables considered inthis study 
include, drum thickness, width, the angle between the hinged 
pin and the end of the Lining, and the distance from the 
pivot to the center of rotation. 

Inequality constraints- Ine-sided conditions which 
must be mathematically satisfied for the design to be 
acceptable. The common mathematical term is G(x <0 or 
G(X) >0. If the inegality condition is satisfied on G(X), the 
d2sign is acceptabl2, (feasible). If it is not satisfied, 
the design is not acceptable (infeasible). Constraints 
considered here incluie, vehicle stopping time, maximum drun 
tamperature, and actuating force. 

Side constraints- Upper and lower bounds on the 
individual design variables xX. The common mathematical 
representation is ion Ki< ae 

Design space- The n-dinensional mathematical space 
Spanned by the vector of design variables X. 

Active constraint- Constraint G. (X) is called active 

if its value is zero (or near oot for computational 


purposes). 
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Phacesve constraint- Constraint 6 (x) is inactive if 
Gmx) <0. 
; Violated constraint- Constraint 6, (X) is violated if 
en (X) >0. 
J 


Ce THE OPTIMIZATION PROCESS 
The general design optimization problem can be stated 
mathematically as follows: Find the set of variables Xo 
meip2....8, Which will 
Minimize F(X ) (35) 
Subject to: 


G(X) < 0 eee a (36) 
xtc x. < x? ia 7 (37) 
a ee 


Vector xX contains the set of independent design 
variables hie a1 een « X may represent, for example 
width, thickness, and angies in the brake Optimization. fThe 
objective function used here is the drum volume. 

Equation (36) defines the inegqality constraints imposed 
on the design. For 2xample, if the temperature on the inner 
drum surface must not exceed a specified value f, he 
associeteli design constraint becomes, in normalized form 

a 


me i < 0 7 
T = 


KB 


The lower and upper bounds on the design variables, 
given by £q. (37) , limit the region over which the functions 
F(X), and G (Xx) are defined. These constraints are often 
referred £5 as Side constraints because they form the sides 
or bounds of the n-dimensional spaced spanned by the design 
Watiables X. 

If all the inequalities of Eqns.(36) and (37) are 
satisfied, the design is said to be feasible; if any of 
these conditions ar2 not satisfied, the design is not 


24 











feasible. If F(X) is a minimum and the design is feasibla, 
it is also optimum, or at least, a relative optimum. Note 
that because the objective and constraints may be nonlinear, 
there may be multiple minima in the design space that cannot 
be identified using current methois. While this is a matter 
Beerconcern, Since it is desired to find the true optinun, 
it must b2 remembered that the same mathematical conditions 
exist if the design process is not automated. However, using 
Optimization techniques, it is a simple matter to restart 
the optimization from several initial points in the design 
Space and thereby improve the probability of obtaining the 
true optinum design, a process that would be quite time- 
consuming in manual design. 

Equations (35)-(37) define the nonlinear constrained 
optimization problem. If Eqs.(36) and (37) are not imposed 
on, the design, the optimization problem is defined by 
Eq.(35) alone and is therefore an unconstrained minimization 
problen. 

Most 2onlinear optimization algorithms update the vector 


of design variables by the iterative relationship; 


meee’ + + (ost (39) 


where q is the iteration number, yYeetor S is the direction 
of search in the design space, and the scalar q is referred 
tO aS a noVe parameter which,together with S$, determines how 
Mach the vector X is changed during the q-th iteration. An 
initial design defined py x must be supplied. The 
Optimization process then proceeds in two steps. FLiEst,; 
the direction S, which improves the design, is found, and 
second, the scalar a, is determined which improves the 
design as much as possible when moving in this direction. 
[The process is repeated until there is no further design 
improvement, indicating that this is the optimum attainable 


25 





design. Por further details see Ref. 7. 


Dee COPES AND SUBROUTINE ANALIZ 

In order +t+o sinplify the use of CONMIN and to further 
aid in the design optimization process a Control Program For 
Eagineering Synthesis, COPES, was developed by Vanderplaats 
fRef. 7]. COPES is the main ovorogram (r2call that CONMIN is 
woeitten in subroutine form). The user must supply an 
analysis subroutine with the name ANALIZ, which will 
calculate the various parameters. This subroutine has three 
segments; INPUT, EXECUTION, OUTPUT. 

All parameters which may be design variables, objective 
functions or constraints are contained in a single labeled 
common block called GLOBCM. 

Copes Terminology 

mime COUPES projram currently provides six specific 
capabilities; 

1. Simple analysis, just as if TOPES was not used. 

2. Optimization-Minimization or maximization of one 
calculated function with limits imposed on other 
functions. 

3. Sensitivity analysis- The effect of changing one or 
more design variables on one or nore calculated 
functions. 

4, Two-variable function space-Analysis for all specified 
combinations of two design variables. 

5. Optinum sensitivity- The same as sensitivity analysis 
except that, at 2ach step, the design is optimized with 
respect to the indapendent dasign variables. 

6. Approximate optimization- Optimization using 
approximation techniques. Usually more efficient than 
Standard optimization for up to 10 design variables or 


if multiple optinizations ar to be performed (Ref. 7}. 
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ny eo ecoen ee rLON OF THE COMPOTER PROGRAM 


A. GENERAL PROGRAN ORGANIZATION 

A functional block diagram o£ the program is presented 
Moe Fig. 4. A general description of the subroutines 
contained in the program is given here. Appendices A 
erzough D discuss the preparation of input data, list the 
laportant computer program nomenclature, and list the 


prograa. 


Be SUBROUTINES 
1. Subroutine ANALIZ 
Subroutine ANALIZ organizes the basic analysis used 
meeene Optimization. It controls the reading of the initial 
design description and calculation of the values of the 
objective function, constraints, and all other parameters 
Necessary to solve the problem. COPFS/CONMIN updates the 
design to Minimize/maximize the objective Lune tron, 
iterating until no further improvement in the objective 
panction is possible without violating one of the 
constraints. COPES/CONMIN calls subroutine ANALIZ to obtain 
the function value iuring the optimization. 
2. Subroutine [INPUT 
[This subroutine reads all input data associated with 
the brake analysis. ENSerEucti ons for problen deck 
preparation are given in appendix B. 
Bee oubroutine TEMPR 
This subroutine calculates the heat transfer 
constants such as the thermal capacity of each node and the 
resistance of each node, determines the time increment in 
order to insure a stabl2a solution, and calculates the rate 
of heat generation. In order to calculate the temperature of 


each node, it calls two subroutin2s. From subroutine BRAK it 
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oSbtains the deceleration needed to calculate the rate of 
haat generated and from subroutine TEMA it obtains the 
temperature ris2 of each node. Then it calculates the 
temperatuc2s during tae time that the brake is not in use. 
This subroutine is also sapable of calculating the 
Mempecatuce rise of a drum when a constant rate of heat 
dissipation is givea. 

4. Subroutine TEMA 

This subroutine calculates th2 temperature of each node. 
As mentioned befor2, tne heat is generated on the inner 
surface, and onthe outer side of the drum the heat is 
dissipated by convection. The formulas used were developed 
by the finite difference method, and are given in section 
mele i 2 

pee Subroutine BRAK 

This subroutine calculates the torque, actuating force, 
and the friction moment of one shoe. it also calculates the 
drum volume and th2 deceleration of the machine. The 
Subroutine takes into consideration a constant friction 
coefficient until a temperature of 90°C is reached anda 
linear decrease in the friction coefficient for higher 
temperatures. More details are given in section II-D-4. 

See subroutine JUIPUT 

This subroutine echos the input data and prints out 

the thermal and mechanical information for the brake. An 
example of the output obtained from this subroutine is shown 


in Table 1 and Table 2. 


28 





Vee teats PROSEEN AND RESULTS 


The computer program was «ested with the data specified 
in Table 1. The objective function which was minimized was 
the volune of the drun material. Design variables were the 
drum width, the angle between th2 hinged pin and the end of 
eae Lining, the batio of the pivot to center of rotation 
distance #9 drum radius, and the drum thickness. The side 


constraints (limits) on the design variables were; 





Xesign Lower Upper 
Variable Bound Bound 
ee (3) Width,b OO 80 mm. 
Aa (5) Theta 2, tee bad. 22 5- Pad. 
Sa (12) Ratio,Rd 0.1 d29 
4. (18) Thickness,tk 4O mm. No bound 


The number in parentheses is the location of the 
variable in the COMMON block in the computer progran. 

Constraints wer2 imposed on the actuating force F, the 
maximum tenperatur3, ee , on the inner surface of the drun 


and stopping time, t. 


Constrained Lower Upper 
Variable Bound Bound 
he oy) Force,F 200.0 N-m 2000.05 N—a 
Pomei2s) Tine, t No bound 12007 «Sec. 
3. (%) Temperature,T No bound 2306.0°6 


The vehicle which weights 25700.0 Newtons is stopped 
four consecutive times froma velocity of 90.0 Km/hr to 
zero, with an acceleration period of 20.0 sec. between 
stops. The values of the design variables and the 


constraints before and after optinization are; 
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Objective 
Punction 
Drum Volume 


Before 


Jptinization 


0.754 E-03 a3 


After 


Ootimization 


0.159 E-02 n° 


Design 

variables 

Width 0.08 a 0.08 o 
Theta 2 2. VOeRad. 1,92 Rad. 
ay = O75 02755 
Thickness 0.010 a 0.020 a 
Constraints 

Actuating 

Force ZE15. 10N 2086.3 N 
Stopping 

time (last stop) 7.904 sec. 7.00 sec. 
Temperature 

after Last stop 348.7°C 229.2 °C 


Note that the objective function increased as a result 

Of optimization. This is because the initial design violated 
constraints on stopping time and naximum temperature. 

Further results are listed in Tables 1 and 2. In 

Mention tO optimization, a sensitivity analysis of the 

design variables and a two-variable function space analysis 

for width and thickness were performed. The graphical 

results are given in Figs. 5 through 17. The results can be 

Summeriz2d as follows; 
a. The effect of changing thea inside drum radius with all 
Other design variables held constant; 

As shown in Figs. 5-7, for small inside drum radii the 

drum temperature is very high. The stopping time is 

long and the torgue is low. Inside drum radii over 130 

mm give reasonable drum temperature and stopping time, 


for the example considered. 
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As seen in Figs. 8-10, the affect of changing the drun 
width with all other design variables held constant is 
the same as described above. 

The effect of changing the drum thickness with all 
other parameters held constant is; For a drum thickness 
up to 6 mm, tne stopping tine and drum temperature are 
considerably high. Over 16 am thickness, the stopping 
tim2 remains almost constant. For a small thickness the 
torque is very low due to the high temperatures. For 
thicknesses over 20 aa, the torgue remains about 
constant. 

The effect of changing the angle between the hinged pin 
and the end of the lining is; For a small 92 angle the 
stopping time is very long because the torgue is low. 
The stopping time becomes r2asonable when 6, >1.8 Rad. 
Qbvidusly there is an increas? in the drum temperature 
as @, increases but the overall change in temperature 
is sm2ll. 

meom che two Variable function space, Fig. 17, it can 
be seen that the constant volume line and the constant 
temperature line ar2 almost parallel, this leads to the 
conclusion that for the cycle taken, the drum is a heat 
Sink, and the amount of heat dissipated by convection 


during this cycle is small. 
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VI. DevemsatURS RALSE =sSmMPLIFIED CALCULATION 


A simplified way of finding the temperature rise of the 


drum is by using the 2quation; 
a 2 zOAT (40) 


and setting Q equal to the amount of heat generated using 
Equation (21) fromesection Sll-0-4. (This equation is in 
common us2 in 2ngineering texts (See, for example, Refs. 2 
and 8). The temperature rise calculated this way, is tha 
average temperature of the drum, and not the temperature on 
the interface, which can be much higher (depending on the 
rate of heat generated). Extreme temperature gradients cause 
distortion and excessive surface wear. Therefore it isn't 
always acceptable to use the simplified formula. From 
experience, it has been found that the surface wear 
increases dramatically as interface temperatures approach 
400 to 509°F (205 to 260 “c), [Ref. 8]. 

A comparison of the temperatures calculated on the inner 
drum surface, outer drum surface, and the average drum 
temperature calculated, using equation (40), is given in 
Fig. 18. The graph shows the temperature rise for a vehicle 
Beopiped from a velocity of 90 kn/hr. From this graph, it can 
be seen that the drum temperatur2, based on equation (40), 
after the vehicle stopped is about the average temperature? 
of the inner and outer surface temperatures. 

The results show that the drum will reach an uniforn 
temperature of about 53°c, in 15 sec, after the vehicle has 
stopped. 

Calculating the tamperature with the simplified formula, 
Can lead t9 errors in the time needed to stop the vehicle. 


Because the temperature calculated with the simplified 
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formula is lower than the temperature at the friction 
M@aterface, the calculated friction coefficient is higher 
than the actual) Gpictron ~ecoetficient. Therefore the 
calculatei stopping time will be shorter than the real 
stopping time. Alb this is trues, Provided the friction 
material behaves as assumed in section I-D-4. 

Because high temperatur? is detrimental to both the 
stopping ability and the wear characteristics of the brake, 
it is important that the interfac2 temperature be calculated 
with reasonable accuracy in design. Fig. 18 clearly shows 
the temperature differences resulting from the two 
approaches. 

This difference in results is compounded when the 
Simplified equation is used for design. Table 3 presents the 
design results based on the simpiified approach. This design 
represents an apparent material savings of 27%. However, 
when this optimum is analysed using the finite difference 
heat transfer solution, the maxinum temperature is 268.8°C 
and the last stopping time is 7.54% sec. This time violates 
mre constraint by 7.7%. Perhaps more importantly, the 
temperature at the interface of about 269°C would surely 
lead to premature failure. Therefore this design is clearly 
too unconservative to be acceptable. 
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Vil. CONCLUSION 


In summary, a numerical optimization program is an 
effective way of finding a solution to an engineering 
problem, provided reasonable care is used in formulating the 


problen. 


VITI. RECOMMENDATIONS FOR FUTURE INVESTIGATIONS 


The study has shown the feasibility of using numerical 
Optimization in the design of Internal-Expanding Rim Brakes 
with two Leading shoes. Further studies on the same design 
may be pursued by eliminating som? of the restrictions. For 
example; 

1. To add heat dissipation by radiation. 

2. To investigate drum temperatures for a drum with fins. 

3. To take into consideration changes in the surface 
pressure as a function of friction coefficient. 

4%, To repeat all the calculations for a drum in which 
there iS one trailing shoe and one leading shoe. 


5. To add the effect of centrifugal forces for clutches. 
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APPENDIX A 
Uist yOr PRRAUETERS 


A complete listing and description of all variables used 
in the ovrogram, is not practical. The variables listed in 
this appendix are common to s2veral subroutines of the 
program and will assist the raader in a study of the 
program. The Global location is the location of the 
parameter in the common block called GLOBCM. This common 
block is the means by which information is transfered 
between the subroutines and the COPES/CONMIN progran. 


Global Fortran Math. Definition 
Location Name Symbol 

1 Ren r Inside drum radius (m) 
2 af c; Torque of one shoe (N-m@) 
3 WDTH b Drum width (nm) 

4 PRSA p Pressure between Lining 


and drum (N/m*) 

=) TETA! 01 The angle between the 
hinged pin and the (Rad.) 
begining of the lining 

6 PETA2 G 2 The angle between the 
hinged pin and the end 
of the lining (Rad.) 


FRMNT Mf Friction moment (N-m) 
ANMRT Mn Normal momemt (N-m) 
ACFRC F Actuating force (N) 

10 C d Distance from actuating 


force to the hinged pin(m) 


G7 





11 
12 


13 
14 
15 


16 
VW 
18 
19 


20 
21 
22 
23 
24 
25 
26 
27 
Zo 32 


3 
34 
55 
36 
a 
38 
Bo 
40 


4] 
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RD 


CAIU 
HM IU 
AMIU 


SREC 
RO 
THK 
DX 


RIIER 
W 

DCCE 
TOT 
ECEN 
NWR 
TIME 
VOL 
TEPL (5) 


NWR 
NWRA 
NWRO 
NEL 
NS EG 
a. 
G 
K 


HCNV 


Sear 


tk 


W 
ge 
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Heat generated (J/sec.) 
Distance Erom pivot to 
center of rotation (an) 
Coma fErmction coectficient 
Howe Ebiceion coefficient 
Fuiction coefficient at 
any temperature 

Drum surface area (m*) 
Outside drum radius (a) 
Drum thickness (m) 

An incremental thickness 
(2) 

Wheel radius (a) 

Car's weight (N) 
Deceleration (m/sec.’ ) 
Total time 
BGentrie ty (mn) 


(sec.) 


Write statement control 
Time (sec.) 

Drum volume (m’) 
Temperature at time prt 
(sec.) 

Write statement control 
Write statement control 
Write statement control 
Number of elements 
Number of segments 
Constane 

Gravitational constant 
Thermal conductivity 
(J/m=C) 

Convection heat coeff. 
(W/m2-C). 

Specific heat (J/Kg-C) 





43 
44 
45 
46-50 
51-57 
pe 6 3 
64 
65 
po- 7 2 
a= 79 
80 
81 
82 


83 
84 
85 
86 


ACOF 
TINT 
ZMAN 
NSHO 


BL 
ro 


ct 
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Density (Kg./m’) 

Time increment (séc.) 
fax. temp. difference (°C) 
Temp. at time p (sec.) 
Heat resistance (°C/J) 
Heat capacity (J/°C) 
Biot moduli 

Fourier moduli 

Control parameter 
Control parameter 
Number of elements+t1 
Number or elsments+2 
The angle at which the 
pressure between the 
lining and drum is 
maximum. (Rad.) 
Constant 

Initial temperature (°c) 
Time increment (sec.) 


Number of shoes 





APPENDIX B 
TNSTRUCTIONS FOR PROBLEM DATA PREPARATION 


Although the procedure 1s straight forward, preparation 
9f input data for the program requires attention. Errors are 
meow tO make and difficult to locate. Input data is 
dzscribed here for the brak2 analysis. For instructions on 
data preparation for optimization see Ref. 7. Input data 
should, in general, follow the steps outlined below. The use 
of the standard FORTRAN Eighty cColumm Coding Sheet is 
recommended. Integer constants nust be right justified in 
the appropriate field. There are eight input cards, read by 
subroutine INPUT, to describe th2 initial design, material 


properties and constants. Card format is given in 
parenthesis followed by specific instructions where 
necessary. 

ieee F2est Card (110) - Duty cycle information. 


Cols 1-10 : Total number of consecutives stops 
and accelerations (NSEG) 
Geese cond Card (110,3F10.0) - Duty cycle information. 
ae Ors 1-10 : Control number. 
1 means-deceleration, 
2 means-brake not in use, 
b. Cols 11-20 : Velocity at start of deceleration. 
GG. Cols 21-30 : The velocity at the end of the 
deceleration. 
d. Cols 31-40 ;: The time the tErake is not in use. 
3. Third Card (5110) - Taarmal analysis information. 
a. Cols 1-10 Number of nodes (NEL). 


b. Cols 11-20 : An integer number that controls the 


amount of printout when detailed 
output is required during the 


vehicle deceleration. The amount of 
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ww @ 


ale 


e. 


CoUSA2 i= 30): 


Coisest-—49 = 


Cols 41-50 : 


lines written, depends on the 
stopping time and time increment. 
(NWR). 

An integer number that controls 

the amount of printout when 
detailed output of the temperatures 
is required during the period that 
the vehicle is not in use. The 
amount of lines written depends on 
the period length that the brakes 
are not used (NWRA). 

An integer number that controls the 
amount of printout when detailed 
output of the temperatures are 
required during the period of 
constant heat generation. (NWRQ). 
Number of braking shoes in the 


machine. 


Fourth Card (7F10.0) - Brake dimensions. 


ae Gors 1-10 &: 

Ds Gots 11-20: 

Se. Cols 21-30 : 

rs Cols 31-40 : 

2. Cols 41-50 3: 

fg Cols 51-60 : 

fem 6 COlS 61-70 ; 

Petth Card (7F10. 
a. Cols 1-10 


Inside drum radius (RI). 

Deum width (WDTH). 

Drum thickness (THK). 

Ratio of distance from pivot to 
center of rotation and inside 
radius (RD). 

Drum density (RHO). 

Angle between hinged pin and the 
beginung of che lining (TETA?) . 
Angle between hinged pin and the 
endsotetne lahing (TETA2). 


0) - Thermal and friction 
ii se wa tion. 
meat conduction coeriicient (K). 


(real number). 


oH 





DeeemcOleelil=20 = Heat convection coefficient (HCNV). 
Gempcols 2l=30) 2 Specizric heat of the drum (SPHT). 
4d. Cols 31-40 : Max. temperature difference between 
Condmerict10n Cosctficient and hot 
friction coefficient (DFT). 
ecolsedd=50 Ss Cold friction coefficient (CMIV). 
mee COlS 51=60 * dot friction coefficient (HMIV). 
g- Cols 61-70 = Initial temperature (TINI). 
6. Sixth Card (2F10.0) - Machine information. 
ae Cols 1-10 : Vehicles weight (W). 
b. Cols 11-20 : Wheel radius (RTIER). 
7. Seventh Card (5F10.0) - Analysis constants. 
a. Cols 1-10 : Maximum pressure between lining 
and drum (PRSA). 
Dawe COlS 11-20 : Constant 3.17415927 
em COlLS 21-30 © Gravitational constant (G). 
d. Cols 31-40 : Increment of time (ZMAN). 
e. Cols 41-50 : The angl2 of maximum pressure 


Seemesloght Card 


Cols 1-10 


(TETAA). 


(CiElO2 0) s= Print control. 


An integer number can be zero or 1. 
I£ zero (or a blank card) - only 
the final results are printed. 

If 1- the temperature at time 


increments are printed. 
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CALL OUTPUT (ICALC) 


APPENDIX D 
PROGRAM LISTING 
SUBROUTINE ORGANIZES THE BRAKE ANALYSIS 
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